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For networks of pulse-coupled oscillators with delayed excitatory coupling, we analyze the firing 
behaviors depending on coupling strength and transmission delay. The parameter space consisting 
of strength and delay is partitioned into two regions. For one region, we derive a low bound of 
interspike intervals, from which three firing properties are obtained. However, this bound and these 
properties would no longer hold for another region. Finally, we show the different synchronization 
behaviors for networks with parameters in the two regions. 
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For decades, complex networks have been focused on 
by scientists from various fields, for instance, sociology, 
biology, chemistry and physics, etc. In particular, net- 
works of pulse-coupled oscillators, as an important class 
of interconnected dynamical systems, have gained in- 
creasing attentions because of their intimate relationship 
to natural systems as diverse as cardiac pacemaker cells, 
flashing fireflies, chirping crickets, biological neural net- 
works, and earthquakes (cf. [11, 0,13]) ■ A pioneering work 
on modeling and analyzing pulse-coupled units was done 
by MiroUo and Strogatz Inspired by Peskin's model 
for self-synchronization of the cardiac pacemaker, they 
proposed a pulse-coupled oscillator model with unde- 
layed excitatory coupling to explain the synchronization 
of huge congregations of South East Asian fireflies. With 
the framework of the MiroUo-Strogatz model, many the- 
oretical and numerical results on pulse- coupled networks 
have been obtained 
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Pulse-coupling is difficult to handle mathematically be- 
cause it introduces discontinuous behavior into the other- 
wise continuous model and so stymies most of the stan- 
dard mathematical techniques 4S]. Particularly for de- 
layed pulse-coupling, the mathematical analysis of collec- 
tive dynamics of networks becomes a challenging prob- 
lem. Past research experience indicates that some under- 
lying facts and assumptions about firing behaviors play a 
crucial role in mathematical analysis [Ijl, 0, H, S 
For example, in 0, synchronization was proved by 
making use of the fact that the firing order of oscillators 
is always preserved for complete and undelayed pulse- 
coupling; in an assumption about firing times made 
the analysis easier by reducing the number of case dis- 
tinctions; in [l^ . the proof of desynchronization was es- 
sentially due to a low bound of interspike intervals. 

In this Letter, networks of pulse-coupled oscillators 



with delayed excitatory coupling are studied. We ana- 
lyze the firing behaviors depending on coupling strength 
and transmission delay. The parameter space consisting 
of strength and delay is partitioned into two regions. For 
one region, we give a low bound of interspike intervals. 
By using the bound, three firing properties are derived, 
which would be very helpful for discussing synchroniza- 
tion of networks and stability of periodic solutions. Un- 
fortunately, these properties no longer hold for another 
region. Furthermore, the different synchronization be- 
haviors for networks with parameters in the two regions 
are presented. 

We consider a system of N identical oscillators which 
are pulse-coupled in a delayed excitatory manner. As in 
Iq^I, the coupling structure is specified by the sets Pre(i) 
of presynaptic oscillators that send pulses to oscillator i, 
or the sets Post(i) of postsynaptic oscillators that receive 
pulses from oscillator i. A phase variable (/'i(t) € [0, 1] is 
used to characterize the state of the oscillator i at time t. 
In the case of no interaction, the dynamics of (jji is given 
by 



d0,(t)/dt = l, 



(1) 



namely, the cycle period of the free oscillator is 1. When 
reaches the threshold (pi = 1, the oscillator i fires and 
jumps back instantly to zero, after which the cycle 

repeats. That is. 



- 1 ^ = 0. 



(2) 
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Because of the transmission delay, the oscillators interact 
by the following form of pulse-coupling: if oscillator i fires 
at time t, it emits a spike instantly; after a delay time r, 
the spike reaches all postsynaptic oscillators j G Post(i) 
and induces a phase jump according to 

0,(i + T) =/-i(min[l,/(0,((t + r)-))+e,,]), (3) 

where Sij > is the coupling strength from oscillator i 
to oscillator j; and the function / is twice continuously 
differentiable, monotonously increasing, /' > 0, concave 
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down, /" < 0, and satisfies /(O) = 0, /(I) = 1. For a 
more detailed introduction of the model, see [1, 0, H B 
[13, In this Letter, we further assume the following: 
(i) The coupled system starts at time t = with a set 
of initial phases < 4>i{Q) < 1; (ii) there is no self- 
interaction, i.e., i ^ Pre(i) for any oscillator i; (iii) < 
r < 1; and (iv) the coupling strengths are normalized 
such that for all oscillator i, J2jePrc(i) ^ji = ^ with < 
e < 1. 

We partition the parameter space A — {{t, e)\Q < t < 
1, < £ < 1} into two regions 

^i = {(T,£)e^|/(r) +£<!}, 
^2 = {(T,£)e^|/(r) +£>!}. 

First of all, we use "proof by contradiction" to prove that 
no oscillator can fire twice in a time window of length t, if 
parameters (r, e) ^ Ai- Let ti and t2 with ti < t2 be two 
successive firing times of oscillator i. Suppose t2 — ti < r. 
We claim that if t2 > t, there must be some oscillator 
i' G Pre(z) which fires more than once in the time interval 
(ti — T,t2 — t] n [0,oo). In fact, this comes from the 
monotony and concavity assumption of the function /. 
Since /' > and /" < 0, we have that for any < S < 1, 
if < Si < 6I2 < - S), then /"H/C^'i) + S) - 0i < 

f^\fi02) + <5) - 02, namely the property {A7) in ^. It 
implies that in the same circle, the later the spike arrives, 
the larger the induced phase jump is 0, Sj- Therefore, if 
all the presynaptic oscillators j € Pre(i) fire at most once 
in {ti — r, t2 — t] n [0, 00), then in the time interval (ti, t2] 
the sum of the phase jumps of oscillator i is not more 
than /~^(/(<2 — ti) +e) — (^2 — ^i), i-e., the sum reaches 
its maximum if all spikes arrive at time ^2 simultaneously. 
It means .^.(ta) < f-\f{t2-h)+e) < r\f{T)+e) < 1, 
which contradicts that oscillator i fires at t2- Thus, there 
exists some oscillator i' E Pre(?) firing more than once in 
(ti-r, i2-T]n[0, 00). Lett3,t4 € (ii - t, t2 - t] n [0, 00) 
with ^3 < t4 be two successive firing times of oscillator i' . 
From t2 — ti < T, it follows that t^ — t^ < t. Similarly as 
above, if t4 > r, then there must be some oscillator i" e 
Pre(i') which fires more than once in the time interval 
(ts — T, f4 — r] n [0, 00). Repeating the derivation leads to 
a finite sequence of pairs of firing times: 

{tl,t2} {t3,t4} > {t2n-l,t2n} (4) 

which satisfies [t2k+i,t2k+2] C (i2fe-i - t, t2k - r] n [0, oo) 
for A: = 1, . . . , n - 1; t2n < T, t2fc > T for fc = 1, . . . , n - 1; 
and each term of ([4]) is two successive firing times of some 
oscillator. Particularly, i2n-i and t2n are two successive 
firing times of some oscillator zq. However, similarly as 
the argument of 4>i{t2) < 1, according to t2n < t, /(t) + 
e < 1 and the assumption that the coupled system starts 
at time i = 0, we can get 4>io{t2n) < 1- It contradicts 
that oscillator io fires at t2n- This contradiction comes 
from our hypothesis t2 — ti < t. For a more detailed 
proof, see [l3|. As a consequence, we get 



Theorem 1: If parameters (t, e) G ^i, all interspike 
intervals of each oscillator in the coupled system must be 
longer than the delay time t. 

Here and throughout, "interspike interval" is referred 
to as the time between two successive firing activities of 
an oscillator. However, as opposed to Theorem 1, at each 
(t, e) G A2, the coupled system has solutions in which 
some interspike intervals are not longer than r. The sim- 
plest example is that the oscillators with initial phases 
0i(O) = ••• = 0jv(O) fire synchronously with a period 
t = r, if (t, e) S A2- In the rest of the Letter, one will 
see that this can cause significantly different dynamical 
behaviors of the system at (r, e) £ Ai and at (r, e) G A2, 
especially the different firing behaviors. Before discussing 
the difference of firing behaviors, let us give some defini- 
tions and notations. Denote the time at which oscil- 
lator i fires its m-th time. Clearly, the firing times f^, 
i ^ 1 . . . , N , m > 1, are determined by initial phases. For 
a given set of initial phases [(^i(0), . . . , 0Ar(O)], the solu- 
tion [01 (i), . . . ,(pN{t)] is said to be a period-d solution 
if there exist a Atg > 0, and positive integers M and d 
such that the firing times of arbitrary oscillator i satisfies 
^m+d ~ = ^'•0 for all TO > M. For a given set of initial 
phases [(/)i(0), . . . , 07v(O)], the solution [(pi{t), . . . ,(j)N{t)] 
is said to be a completely synchronized solution if there 
exists a T > such that the phase variables of arbitrary 
oscillators i and j satisfy tpiit) — fj{t) for all t >T. For 
the convenience of later use, we let Sij — for j ^ Post(i). 
Then, the phase jump ([3|) also holds for j ^ Post(z). 

By using Theorem 1, we conclude that if (r, e) e Ai, 
any solution of the coupled system possesses the following 
properties: 

Property 1: For oscillators i and j satisfying Sij — eji 
and Eki = Ekj for aU k e Kij := {1, 2, . . . , TV} \ {i,j}, if 
C < tL,i then tj„^+i < t^^.+i, i.e., the firing order of i 
and j is always preserved. 

Property 2: For oscillators i and j satisfying Sij — Sji 
and Ekt = Skj for all k e K^j, if — P^^, then ^^(i) = 
(f>j{t) for all i > C- 

Property 3: If [(j)i{t), . . . ,(f)M[t)\ is a completely syn- 
chronized solution, then it is a period-one solution. 

In fact. Theorem 1 implies that in the case of (t, e) G 
Ai, a spike of oscillator i must reach oscillators j G 
Post(i) before i emits the next spike. Thus, for oscillators 
i and j satisfying = Sji and eki = £kj for all k G i^ij, 
the instantaneous synchronization i^. = tl-^_ can lead to 
(j)i{t) = <t>j{t) for all t > t^j. (Property 2). For the same 
reason, if [0i(t), . . . ,4>Nit)] is a completely synchronized 
solution, then tl^_^i —tin = 1~ [/~"^(/(''') + ^)~'''] lor all i 
and tin — T ■ That is to say, any completely synchronized 
solution is a period-one solution with the final interspike 
interval being 1 — [/"^(/(r) + e) — t] (Property 3). 

Due to space limitations, here we prove Property 1 for 
the case of iV = 2. By Property 2, we only need to prove 
that ii t]n^ < i^^, then t^^.+i < C^+i- The proof is 
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divided into four cases: 
Case 1: > C+i- 
In this case, we have fi. 



(a) 



-mi + i — ''ma ^ ''ma + l- 

Case 2: t^, + r < t^^ 

In this case, since (i^a' 
oscillator 2 cannot receive any spikes from oscillator 1 
in the time interval (^maj^mi+i]- This, combined with 
,)+) < Mtl^J < 1, leads to Mt) < Mt) for 



^ '-mi + l- 



2 

ma > 



O = 02((t 

ah t e 
Case 3 



'ma ' 



mi+1 



It implies i^^^+i < tl^^+i- 



Since (t^ 



mi ■ T and 01 (i^ 



'/'2(ii 



mi ' 

r) < 



' : '' m 
'^1(4 



c 



(^mi 



''"i^mi + 1 + 



and 



t), similarly as Case 2 we can 
et (t)2{t) < 01 (t) for all t e (ii„^ + T,tmi+i]- It implies 



'•mi + l ^ ''ma + 1- 

Case 4: < 

^2 ,ti ^ 



<02(4i +^)- 
, by Theorem 1 



V'ma ' mi 
'/'l(i™a) < 



T and 01 (t^^ 
Since (CC, +t) C (4i,4i - 
oscillator 2 cannot receive any spikes from oscillator 1 in 
2 '1 -Lr). This, combined with 02((ima)^) < 
1, leads to 02((4i + ^)") < '^i((4i + ^)")- 
Let/o = /(0i((4i+T)-))-/(02((4i+r)-)). Because 
the spike emitted by oscillator 1 at reaches oscillator 
2 at + T, we have /(02(tLi + r)) - /(0i(4i +^)) = 
£i2 — /o- It can be claimed that /(02((ima +''')")) ~ 
/(0i((t^2 +t)~)) < £12 — /o- Indeed, this comes from the 
property (A5) in i: /(02)-/(0i) > /(^2+^)-/(ei+^), 
if 6*1 < 6*2 and < (5 < 1 - 6*2 . Denoting At = 
from the property (A5) in Q we get /(02((ima +''')")) 

/(<^l((4a+^)")) - /(02(Ci+T) + At)-/(0l(4i+T) + 

At) < /(02(4i+r))-/(0i(C^+r)) = £i2-/o. Because 
the spike emitted by oscillator 2 at reaches oscillator 
1 at ^2^^ +T, we have /(0i(t^, +t)) = min[l, /(0i((t^, + 
t)^))+£2i]. So, if /(0i(t^2+T)) = 1 , then from Theorem 

1 it follows that /(02(t^a + ^)) < 1 = /(0l(Ca + ^)); 

/(</'i(^ma + ''')) < I' then from the above claim it follows 

that/(02(4a+^))-/('^l(*ma+^)) = /('/'2((4a +^)"))- 
/(</>! ((ima + ^)")) - £21 < £l2 " /o - ^21 = -/o < 0. It 

implies ti^^+i < C+i- 

In fact, we proved that for the case of = 2, if t\ < 



^mi ' 



^ma' then tmj^_i 



< t: 



ma + l- 



For the case of > 2, the 
proof is similar, and also can be divided into the above 
four cases. The distinction is that when N > 2, t^i+i = 
may happen in Cases 2-4. This derives from the 
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ma + l 

fact that two oscillators are likely to desynchronize, while 
the other oscillators try to synchronize them 

Numerical analysis shows that from any initial phases, 
the coupled system approaches a period solution with 
groups of synchronized oscillators 0, i S El, El. In 



larger networks, the oscillators can be divided into groups 
in a combinatorial number of ways, and exponentially 
many periodic solutions are present 0], which greatly 
increases the complexity of firing behaviors. Proper- 
ties 1-3 indicate that the firing behaviors of the cou- 
pled system at (r, £) £ Ai are relatively simple. How- 
ever, when parameters (r, £) G A2, there may be some 
solutions, which do not possess some or all of Proper- 
ties 1-3. It makes firing behaviors more complicated. 
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FIG. 1: Firing times of four all-to-all pulse-coupled oscilla- 
tors with T = 0.9 and e = 0.6. Vertical dashed lines are 
used to indicate the boundaries of periods, (a) Initial phases 
[0i(O), 02(0), 03(0), 04(0)] = [0.1766,0.4298,0.4079,0.7061]. 
Two periods (i) and (ii) are presented, (b) Initial phases 
[0i(O), 02(0), 03(0), 04(0)] = [0.4974,0.2492,0.8932,0.8501]. 
Four periods (i), (ii), (iii) and (iv) are presented. 



Whether or not such solutions exist depends on param- 
eters (t, e) and coupling strengths £,jj . For the system 
with (t,£) e A'^ := {(r,£) £ -4|/(r) + e > 1} and 
Eij — e/{N—l), i ^ j (hereinafter referred to as all-to-all 
coupling), such solutions always exist. Furthermore, for 
any such solution, there must be some interspike intervals 
not exceeding the delay time r. Otherwise, by previous 
arguments, the solution possesses Properties 1-3. As an 
example, we simulate a network of A^ = 4 all-to-all cou- 
pled oscillators with t = 0.9, e = 0.6. We use for / an 
example of the leaky integrate-and-fire model 
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where / = 1.05. In Fig. [TJa), a period-four solution is 
given. In this solution, the firing order of oscillators 1,2 
(or 3, 4) is not preserved, e.g., < i| but > t%\ and the 
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does not mean 
tl but t\>tl. In 
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instantaneous synchronization t\- 
(j>i{t) = 0j(t) for all t>t\^., e.g. 
Fig. [ijb), a period- two completely synchronized solution 
is given. In addition, one can see that in Fig. [TJa) and 
(b) , most interspike intervals of the oscillators are shorter 
than the delay time r = 0.9. 

Completely synchronized solutions, as a special type 
of periodic solutions, have been widely studied [1, ^, 
0, H, [13, [m. Ell- The following analysis demonstrates 
the different synchronization behaviors for networks with 
(t, e) G ^1 and (t, £) S ^2- In [l2|, we proved that 
under the assumption /(2r) -f £ < 1, from any initial 
phases (other than 0i(O) = • • • = 0jv(O)), all-to-all pulse- 
coupled oscillators with delayed excitatory coupling can- 
not achieve complete synchronization. In factj_we can 
extend this result to the case of (r,£) e Ai (see [ij). In- 
terestingly, we found that when parameters (r, e) S A2^ 
completely synchronized solutions become prevalent. In 
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FIG. 2; Dependence of PN{T,e) on network size A'^. (a) r — 
0.55, e = 0.4. (b) r = 0.7, e = 0.35. (c) r = 0.8, e = 0.3. 



reported. For the species P. malaccae, the normahzed de- 
lay is « 0.36, and perfect synchrony is usually achieved. 

In summary, our analysis demonstrates different dy- 
namics for pulse-coupled networks with (r, e) G Ai and 
(t, e) G Ai. For the region Ai, we derive a low bound 
of interspikc intervals and three firing properties, which 
provide a basis for future researches addressing the dy- 
namics in networks, e.g., stability of periodic solutions. 
The difference of synchronization presented at the end 
of the Letter is useful for understanding and interpreting 
synchronization phenomena in some natural systems. 

This work was supported by the National Sci- 
ence Foundation of China under Grants 60574044 and 
60774074. 
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In Lemma 1 of [l^j, we prove that if the coupling is all- 
to-all and parameters r, e satisfy /(2r) -I- e < 1, then 
all interspike intervals of each oscillator in the coupled 
system are longer than 2r. For normalized coupling and 
parameters (r, e) € Ai, the proof is simalar. 
The proof is almost the same as that in [T2| |. Here, we 
briefiy describe the proof process. By Theorem 1 in this 
Letter, Lemma 1 in [l2| becomes that if an oscillator 
fires at time t\ and t2 with t\ 7^ t2, then \ti —t2\ > t. Al- 
though the result of Lemma 1 is weakened, Lemmas 2-7 
and Theorem 1 in fl^ still hold. Moreover, all the deriva- 
tions need not be changed except that of Lemma 7. For 
the proof of Lemma 7, an additional but straightforward 
case distinction is required. 



